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Abstract. We consider the free boundary problem for the plasma- vacuum interface in ideal compress- 
ible magnetohydrodynamics (MHD). In the plasma region the flow is governed by the usual compressible 
MHD equations, while in the vacuum region we consider the pre-Maxwell dynamics for the magnetic 
field. At the free-interface we assume that the total pressure is continuous and that the magnetic field 
is tangent to the boundary. The plasma density does not go to zero continuously at the interface, but 
has a jump, meaning that it is bounded away from zero in the plasma region and it is identically zero 
in the vacuum region. Under a suitable stability condition satisfied at each point of the plasma-vacuum 
interface, we prove the well-posedness of the linearized problem in conormal Sobolev spaces. 



1. Introduction 
Consider the equations of ideal compressible MHD: 

dtp + div (pv) = 0, 

dt{pv) + div {pv ® V - H H) + Vq = 0, 
dtH X {vxH) = 0, 
L dt{pe+\\H\'^)+div{{pe+p)v + Hx{vxH)) = 0, 



(1) 



where p denotes density, u G R'^ plasma velocity, H ^ magnetic field, p = p{p, S) pressure, q = 
p + total pressure, S entropy, e ^ E + total energy, and E = E{p, S) internal energy. With 

a state equation of gas, p = p{p, 5), and the first principle of thermodynamics, (f ) is a closed system. 
System (1) is supplemented by the divergence constraint 



div H = 



(2) 



on the initial data. As is known, taking into account (2), we can easily symmetrize system (f ) by rewriting 
it in the nonconservative form 



p dt dt 

^-{H ■V)v + Hdivv = 0, ^=0, 
dt dt 



(3) 
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where pp = dp/ dp and A/dt = 9^ + (w • V). A different symmetrization is obtained if we consider q instead 
of p. In terms of q the equation for the pressure in (3) takes the form 



dt 



dH 
lit 



div f = 0, 



(4) 



where it is understood that now p ~ p{q — p/2, S) and similarly for pp. Then we derive drvv from (4) 
and rewrite the equation for the magnetic field in (3) as 

dH , , p„ 
- [H ■ \J)v - 



dt 



P 



d,_^ d^l 
dt J 



dt 



0. 



(5) 



Substituting (4), (5) in (3) then gives the following symmetric system 



Pp/P 

iPp/p)H'' 






ph 
O3 




-iPp/p)H 
O3 

{pp/p)H (2 




H 



0^ 
1/ 



V 



ipp/p)v 

V 

iPp/p)H^v 




V 



V 

H 



V 



(6) 



V- ~{pp/p)Hv 
pv ■ V/3 ~H ■ V/3 

-H-Vh {h + {Pp/p)H®H)v^ 






0^ 









V 




H 




\sj 



where = (0, 0, 0). Given this symmetrization, as the unknown we can choose the vector U = U{t, x) = 
{q,v, H, S). For the sake of brevity we write system (6) in the form 

3 

Ao{U)dtU + Y,Aj{U)djU = Q, (7) 

i=i 

which is symmetric hyperbolic provided the hyperbolicity condition Aq > holds: 

P>0, pp>0. (8) 

Plasma- vacuum interface problems for system (1) appear in the mathematical modeling of plasma 
confinement by magnetic fields (see, e.g., [9]). In this model the plasma is confined inside a perfectly 
conducting rigid wall and isolated from it by a vacuum region, due to the effect of strong magnetic fields. 
This subject is very popular since the 1950-70's, but most of theoretical studies are devoted to finding 
stability criteria of equilibrium states. The typical work in this direction is the classical paper of Bernstein 
et al. [2]. In astrophysics, the plasma- vacuum interface problem can be used for modeling the motion of 
a star or the solar corona when magnetic fields are taken into account. 

According to our knowledge there are still no well-posedness results for full [non- stationary) plasma- 
vacuum models. More precisely, an energy a priori estimate in Sobolev spaces for the linearization of a 
plasma- vacuum interface problem (see its description just below) was proved in [21], but the existence of 
solutions to this problem remained open. In fact, the proof of existence of solutions is the main goal of 
the present paper. 

Let 17+ (t) and fl~{t) be space-time domains occupied by the plasma and the vacuum respectively. 
That is. in the domain f2"'"(t) we consider system (1) (or (7)) governing the motion of an ideal plasma 
and in the domain fl^ (t) we have the elliptic (div-curl) system 

VxH = 0, div-H = 0, (9) 

describing the vacuum magnetic field TL G K'^. Here, as in [2. 9], we consider so-called pre-Maxwell 
dynamics. That is, as usual in nonrelativistic MHD, we neglect the displacement current (l/c)dtE, 
where c is the speed of light and E is the electric field. 

Let us assume that the interface between plasma and vacuum is given by a hypersurface T{t) = 
{F{t,x) ~ 0}. It is to be determined and moves with the velocity of plasma particles at the boundary: 



dF 

-^0 onr(0 



(10) 
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(for all < G [0,T]). As _F is an unknown of the problem, this is a free-boundary problem. The plasma 
variable U is connected with the vacuum magnetic field Ti through the relations [2, 9] 

[q]=0, H-N = 0, n-N^O, onr(t), (11) 

where N — VF and [q] = q\r — ^I'Hj^p denotes the jump of the total pressure across the interface. These 
relations together with (10) are the boundary conditions at the interface T{t). 

As in [11, 20], we will assume that for problem (1), (9)-(ll) the hyperbolicity conditions (8) are 
assumed to be satisfied in 51+ (t) up to the boundary r{t), i.e., the plasma density does not go to zero 
continuously, but has a jump (clearly in the vacuum region n~{t) the density is identically zero). This 
assumption is compatible with the continuity of the total pressure in (11). 

Since the interface moves with the velocity of plasma particles at the boundary, by introducing the 
Lagrangian coordinates one can reduce the original problem to that in a fixed domain. This approach 
has been recently employed with success in a series of papers on the Euler equations in vacuum, see 
[5, 6, 7, 8, 11]. However, as, for example, for contact discontinuities in various models of fiuid dynamics 
(e.g., for current-vortex sheets [3, 19]), this approach seems hardly applicable for problem (1), (9)-(ll). 
Therefore, we will work in the Eulerian coordinates and for technical simplicity we will assume that the 
space-time domains D,^{t) have the following form. 

Let us assume that the moving interface r(t) takes the form 

r(<) ~ {{xi,x') e XI = ^{t,x')}, 

where t e [0,T] and x' = (x2,X3). Then we have r2^(t) = {xi ^ (p{t,x')}. With our parametrization of 
r(i), an equivalent formulation of the boundary conditions (10), (11) at the interface is 

dt^^VN, [q] = 0, Hn^O, nN=0 onr(i), (12) 

where = v ■ N , Hn = H ■ N , -Rn = Ti ■ N , N = {1, ~d2f, -d^ip). 
System (7), (9), (12) is supplemented with initial conditions 

[/(o, x) ^Uo{x), xe n+{o), ^(0, x) = Mx), xe r, 

n{o,x) = no{x), X e n-{o), ^ ^ 

From the mathematical point of view, a natural wish is to find conditions on the initial data providing 
the existence and uniqueness on some time interval [0, T] of a solution (11,1-1, ip) to problem (7), (9), (12), 
(13) in Sobolev spaces. Since (1) is a system of hyperbolic conservation laws that can produce shock 
waves and other types of strong discontinuities (e.g., current-vortex sheets [19]), it is natural to expect 
to obtain only local-in-timc existence theorems. 

We must regard the boundary conditions on H in (12) as the restriction on the initial data (13). More 
precisely, we can prove that a solution of (7), (12) (if it exists for all t <E [0,T]) satisfies 

divi? = mn+{t) and Hn = on r(i), 

for all t G [0,T], if the latter were satisfied at t = 0, i.e., for the initial data (13). In particular, the 
fulfillment of div H = implies that systems (1) and (7) are equivalent on solutions of problem (7), (12), 
(13). 

1.1. An equivalent formulation in the fixed domain. Let us denote 

:= n {xi ^ 0} , r := n {xi = 0} . 

We want to reduce the free boundary problem (7), (9), (12), (13) to the fixed domains For this 
purpose we introduce a suitable change of variables that is inspired by Lannes [10]. In all what follows, 
H'^{ijj) denotes the Sobolev space of order s on a domain oj. The following lemma shows how to lift 
functions from F to M.^ . An important point is the regularization of one half derivative of the lifting 
function w.r.t. the given function (p. For instance, there is no such regularization in the lifting function 
chosen in [12, 13]. 



4 



P. SECCHI AND Y. TRAKHININ 



Lemma 1. Let m > 3. For all e > there exists a continuous linear map ip G H™ ''•^(M^) h- > 4* £ 
ff"(]R3) such that *(0,a;') ip{x'), di'i{0,x') = on T, and 

||9i^'1|l~(r3) < e\\|p\\H^R^)■ (14) 

We give the proof of Lemma 1 in Section 10 at the end of this article. The fohowing lemma gives the 
time-dependent version of Lemma 1. 

Lemma 2. Let m > 3 be an integer and let T > 0. For all e > there exists a continuous linear map 
ip e n^-o^C^([0,r];i7"-J-0-5(R2)) ^ e n]lro^C^i[0,T];H"'-i(R^)) such that *(t,0,x') = (/?(t,a;'), 
di'i/{t,0,x') = on F, and 

!|5l^'|lc([0,T];L = (K3)) < e ||</'llc([0,T];_ff2(R2)). (15) 

Furthermore, there exists a constant C > that is independent of T and only depends on m, such that 

v<^en™-oic^([o,r];iJ"-J-o-5(R2))^ Vj = o,...,™- l, Vie[0,T], 

\\dl^{t, •)|U™-.(R3) < C \\dMt, •)||h..-.-0.5(r2) . 

The proof of Lemma 2 is also postponed to Section 10. The diffcomorphism that reduces the free boundary 
problem (7), (12), (13) to the fixed domains is given in the following lemma. 

Lemma 3. Letm>3 be an integer. For all T > 0, and for all ip G n^1^^C-'([0, T]]H"'-'}-^-^{M?)), satis- 
fying without loss of generality \\ip\\c{[o.T]:H2{R'^)) < 1? there exists a function \1/ & n™LQ^C^([0, T]; 7?™~-'(R'^)) 
such that the function 

<i>{t,x) -.^ {xi + ^{t,x),x') , (i,x) e [0,r] X R3 , (16) 

defines an H"' -diffcomorphism of for all t e [0,T]. Moreover, there holds d{{^- Id) e C{[Q,T]; H"'~^ (M.^)) 
for j = 0,...,m- I, $(t,0,a;') = {p{t,x'),x'), di^{t,0,x') ^ (1,0,0). 

Proof of Lemma 3. The proof follows directly from Lemma 2 because 

di^i{t,x) = 1 + ai*(<,a;) > 1 - ||ai*(t, •)||c([o,T];L~(R3)) > 1 - (-\W\\c([Q.T],m(K^)) > 1/2, 

provided e is taken sufficiently small, e.g. e < 1/2. The other properties of $ follow directly from Lemma 
2. □ 

We introduce the change of independent variables defined by (16) by setting 

U{t, x) U{t, $(t, x)), n{t, x) := U{t, $(t, x)). 

Dropping for convenience tildes in U and %, problem (7), (9) (12), (13) can be reformulated on the fixed 
reference domains as 

P(t/,*)=0 in[0,T]xr2+, V(-H,*)=0 in[0,T]xO-, (17) 

B(C/,-H,^) on[0,r]xr, (18) 

(f/,H)|t=o = (t/o,Ho) inr!+xr!-, v|t=o = '^o on F, (19) 
where P(f/, *) = P([/,*)J7, 

P{U, *) = Ao{U)dt + Ii ([/, *)5i + A2 (U)d2 + A3iU)d3, 
A,iU,^) = -— (a, (U) - Ao iU)dt^ - E iU)dk^) , 



VCH,^') = 



k=2 

V X i5 ' 
div f) 

= (Hiai$,H,,,H.3), f) = (HAr,H25i$i,H3ai$i), 
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B([/,H,^)=( [q] |, \q]=q\..,^^-^\H\l^ 



=0i 



'^N\xi=0 

= vi — V2d2'i' — wa^a^'. 

To avoid an overload of notation we have denoted by the same symbols wat, Hat here above and vm,'Hn 
as in (12). Notice that i'iv|a;i=o = vi~V2d2f~v^d3Lp, T~l-N\xi=o = 'Hi — 'H292</' — 'Hs^av?, as in the previous 
definition in (12). 

We did not include in problem (17)-(19) the equation 

div/i=:0 in[0,r]xl7+, (20) 

and the boundary condition 

Hn = on[0,T]xr, (21) 

where h = (Hn , H2di^i, H^di^i), Hm — Hi ~ iJ2(?2^' — H^d^'^, because they are just restrictions on 
the initial data (19). More precisely, referring to [19] for the proof, we have the following proposition. 

Proposition 4. Let the initial data (19) satisfy (20) and (21) for t = 0. If {U,H,(p) is a solution of 
problem (17)-(19), then this solution satisfies (20) and (21) for all t £ [0,r]. 

Note that Proposition 4 stays valid if in (17) we replace system P{U,'^) = by system (1) in the 
straightened variables. This means that these systems are equivalent on solutions of our plasma-vacuum 
interface problem and we may justifiably replace the conservation laws (1) by their nonconservative form 
(7). 

2. The linearized problem 

2.1. Basic state. Let us denote 

:= {-oo,T]x n^, WT := (-oo,r] X r. 

Let 

{U{t,x),n{t,x),f>{t,x')) (22) 

be a given sufficiently smooth vector-function with U ~ {q, v, H, 5), respectively defined on Qj, Q^, ujt, 
with 

(23) 

ll'^llc([0,T];ff2(K2)) < 1, 

where A" > is a constant. Corresponding to the given (p we construct 4* and the diffeomorphism $ as 
in Lemmata 2 and 3 such that 

> 1/2. 

We assume that the basic state (22) satisfies (for some positive pa, pi G K) 

p{p, S)>po> 0, pp{p, 5) > pi > in Q+, (24) 

dtH + \{w -W)!! - {h-\/)v + HdWu] =0 in Q+, (25) 

VxJ^ = 0, divi) = ing:^, (26) 

dt(p — VN = Q, Hn = on wt, (27) 
where all the "hat" values are determined like corresponding values for {U,H, f), i.e. 

= (nidi^iSr.Srs), ^ = (Hjv,H29i$i,H3ai$i), h = {H N , H2dl^U HaOi^l) , 

p = q-\H\^/2, WAT ui - -Oa^s* - wa^a^', TiN = A - ^■282'^ - ^.383^ , 

and where 

U = {vN,V2dl^l,V3di^i), w ^u - (St*, 0,0). 
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Note that (23) yields 

|lVt,a;*|lvF2,oc([o,T]xR3) < C{K), 

where Vt,x = {dt, V) and C = C{K) > is a constant depending on K. 
It follows from (25) that the constraints 

div/i = inQ+, ^jv = on wt, (28) 

are satisfied for the basic state (22) if they hold at i = (see [19] for the proof). Thus, for the basic state 
we also require the fulfillment of conditions (28) at t = 0. 

2.2. Linearized problem. The linearized equations for (17), (18) read: 

r(C/,$)(<5C/,(5*) := ^nUe.-^e)\e=o = ./ m Qj, 
de 



V'CH,*)((5-H,(5*) 



de 



on ujt, 



where = U + e SU, T-L^ = % + e SJi, ip^ — (p + eSip; 5'^ is constructed from S^p as in Lemma 2 and 
= * + e5^. 

Here we introduce the source terms / = (/i, • . • , /s), ^' = (X;S), x = (xi7X2,X3), and g = (31,32,53) 
to make the interior equations and the boundary conditions inhomogeneous. 
We compute the exact form of the linearized equations (below we drop 5): 



^'{U, $)(t/, *) = P{U, $)f/ + C(C/, $)J7 - {L{U, $)*} 

( 

V'CH,$)(H,«') = VCH,$) + 



V* 



= g\ 



I 



( 

V X -H3 

V V ^2 

q-H-H 

V Hn - H2d2f - Hadstp 
where q := p + H ■ H , vn := vi — V2d2'if — v^ds"^, and the matrix C{U, ^') is determined as follows: 

c(f/, $)r = (r, VyAoiu))dtU + (r, VyMu, ^))diU 

+ {Y,yyA2{U))d2U + {Y,\/yA3{U))d3U, 

8 



l'{U,H,p){U,H,p) 



|a:i=0 



{Y,VyAiU)) :=^2/,; 



fdA(Y) 



Y=U 



Y = {yi,...,ys). 



Since the differential operators F'{U, ^I') and V{H, 4*) arc first-order operators in ^, as in [1] the linearized 
problem is rewritten in terms of the "good unknown" 



U ■.= u 



^1 



diU, H:=H 



^1 



diH. 



(29) 



Taking into account assumptions (27) and (26) and omitting detailed calculations, we rewrite our lin 
earized equations in terms of the new unknowns (29): 



PiU, $)f7 + C{U, - 9i{L([/, $)}=/, 
Y{H,^) = g'. 



(30) 
(31) 
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i'{u, n, 0){u,n, ^) M'{u, n, 0){u, n, ^) 

dt^p + V2d2<p + v^d^ip - vn - vdiVN 
q - H ■ n + [diq](p 
"Hat - d2{'H2^) - d3{'H3<p) 

where vn = vi ~ ^292^ — vads"^, T-Ln = Hi — 'H2d2^ — Ha^a^', and 

[diq] = (5i<7)Ui=o - {H ■ 5iH)U,=o. 

Wc used the last equation in (26) taken at xi = while writing down the last boundary condition in 
(32). ^ 

As in [1, 4, 19], we drop the zeroth-order term in ^' in (30) and consider the effective linear operators 

P^(C/, $){/ := P(C/, $){/ + C(C/, $)[> = /. 

In the future nonlinear analysis the dropped term in (30) should be considered as an error term. The 
new form of our linearized problem for ([/, ip) reads: 

3 

Aodtlf + ^AjdjU + CU = f inQ+, (33a) 

V X X, div f) = S in Q;^, (33b) 

dtip = WW - V2d2ip - hds^ + tpdiVN + gi, (33c) 

q = n-n-[diq]ip + g2, (33d) 

"Hw = 92(1^2^) + dsiHs-p) + 93 on ujt, (33e) 

{U,n,Lp)=0 fort<0, (33f) 

where 

Aa=:Aa{U), a ==0,2, 3, li =: Ii([/, $), C;=C(C/,$), 

h = {Uldi^uUr.Mr,), \) = (?iw,H2C'i$i,H35l$l), 

Hn ^ Til- n2d2^ -Tisds^, Hr^=Hid^^ +m, i = 2,3. 
The source term x of the first equation in (33b) should satisfy the constraint div x = 0. For the resolution 
of the elliptic problem (33b), (33e) the data 2,33 must satisfy the necessary compatibility condition 

Edx^ / gsdx', (34) 
o- Jr 

which follows from the double integration by parts 

/ Edx^ [ dwi)dx= [ fii dx' = [ {82 {H2^) + d3 {ils^) + 53} dx' = [ 53 dx' . 
Jn- Jn- Jr Jr Jr 

We assume that the source terms /, x, S and the boundary datum g vanish in the past and consider the 

case of zero initial data. We postpone the case of nonzero initial data to the nonlinear analysis (see e.g. 

[4, 19]). 

2.3. Reduction to homogeneous constraints in the "vacuum part". We decompose H in (33) as 
v. = H' + v." (and accordingly ^ S)' + S)", f) = f)' + [)"), where V." is required to solve for each t the 
elliptic problem 

Vxi-y^x, div()" = S in^!-, 

VI ^n'l^^gs on r. ^^^^ 

The source term x of the first equation should satisfy the constraint div x = 0. For the resolution of (35) 
the data S, 53 must satisfy the necessary compatibility condition (34). By classical results of the elliptic 
theory we have the following result. 
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Lemma 5. Assume that the data (x, ^,(73) in (35), vanishing in appropriate way as x goes to infinity, 
satisfy the constraint div^ — and the compatibility condition (34). Then there exists a unique solution 
%" of (35) vanishing at infinity. 

Remark 6. In the statement of the lemma above we intentionally leave unspecified the description of 
the regularity and the behavior at infinity of the data and consequently of the solution. This point will be 
faced in the forthcoming paper on the resolution of the nonlinear problem. 

Given H", now we look for T-L' such that 

Vxro'==0, divf)'=0 inQ-, 

q = n-n' -[d^q]ip + g}, (36) 
n'j^ = d2{'H2f) + d'iiUs.if) on cjt, 

where we have denoted <?2 = 52 + • "H". If Ti." solves (35) and H' is a solution of (36) then T-l = H' + H" 
clearly solves (33b), (33d), (33c). 

From (33), (36), the new form of the reduced linearized problem with unknowns (17,1-1') reads (we 
drop for convenience the ' in %' ,g2) 

3 

Aodtif + ^AjdjU + CU ^ f ing+, (37a) 

V X .f) = 0, div [) = in Q^, (37b) 

dttp = vn - V2d2<p - hds^p + tpdiVN + gi, (37c) 

q^n-n-[diq]ip + g2, (37d) 

Hn = d2{'H2'f) + 93(^3'/') on ujt, (37e) 

(C/,-H,¥5)=0 fori<0. (37f) 

2.4. Reduction to homogeneous constraints in the "plasma part". From problem (37) we can 
deduce nonhomogcncous equations associated with the divergence constraint div h = and the "redun- 
dant" boundary conditions Hj\j\^-^^q — for the nonlinear problem. More precisely, with reference to [19, 
Proposition 2] for the proof, we have the following. 

Proposition 7 ([19]). Let the basic state (22) satisfies assumptions (23)-(28). Then solutions of problem 
(37) satisfy 

div h = r inQ^, (38) 
H2d2^ + H^d^ip — IIn — ^diH^ = R on lot. (39) 

Here 

h = {HN,H2di^i,H3di^i), Hn = Hi- H2d2^ - Hsds^- 

The functions r = r{t,x) and R ~ R{t,x'), which vanish in the past, are determined by the source terms 
and the basic state as solutions to the linear inhomogeneous equations 

9ta H — {w -^10 + aAbfu} = Fh inQi,, (40) 

ai$i 

dtR + V2d2R + V3d3R+ id2V2 + dsvs) R= Q onujT, (41) 
where a = r/di^i, Fh = (div //f)/(9i$i, 

fH = ifNjeJr), fN^f5~ f6d2^ - fjds^, Q = {d2{H2gi) + dsiHsgi) - /jvI^^q. 
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Let us reduce (37) to a problem with homogeneous boundary conditions (37c), (37d) (i.e. gi = .92 = 0) 
and homogeneous constraints (38) and (39) (i.e. r = R = 0). More precisely, we describe a "lifting" 
function as follows: 

U = (g,wi,0,0,H,0), 

where q — 32, vi — — gi on lot, and where H solves the equation for H contained in (37a) with v = 0: 

dtH + {{w -V)!! - {h-\/)v + Hdivw] ^ Jh in Qt, (42) 

where h ^ {Hi - H2d2^ - H^d^^, H2,H3), fu = (/s, /e, h)- It is very important that, in view of (27), 
we have wi|a;i=o = 0; therefore the linear equation (42) does not need any boundary condition. Then the 
new unknown 

U'^=U-IJ, H^=n (43) 
satisfies problem (37) with f = F, where 

F^{Fi,...,Fs) = f-V'^{U,^)U. 

In view of (42), Fh = (^5,^5,^7) = 0, and it follows from Proposition 7 that satisfies (38) and (39) 
with r = R = 0. 

Dropping for convenience the indices in (43), the new form of our reduced linearized problem now 
reads 

3 

AodtU + ^AjdjU + CU = F in 0+ , (44a) 

V X i3 = 0, div [) = in , (44b) 

dtip = vn - V2d2ip - vsdsip + ipdiiiN, (44c) 
q^n-n-[diq]ip, (44d) 

'Hn = d2{H2^p) + d3{H3(p) on ujt, (44e) 

iU,n,tp)^0 fort<0. (44f) 

and solutions should satisfy 

div/i = inQ+, (45) 

Hr^ = H2d2ip + H^d^ip — ip OiHn ohojt- (46) 
All the notations here for U and H (e.g., h, Sj, (), etc.) are analogous to the corresponding ones for U 
and H introduced above. 

2.5. An equivalent formulation of (44). In the following analysis it is convenient to make use of 
different "plasma" variables and an equivalent form of equations (44a). We define the matrix 



It follows that 



,0 



-52$ 


-93$ 












h 


= (Hn, 



u {vN, V2di^i,v3di^i) ^fjv, {HN,H2di<S>i,H3di^i) ^fiH. (47) 

Multiplying (44a) on the left side by the matrix 



/I \ 

0^ 77 O3 0^ 

0^ O3 fi 0^ 

V o"^ 0^ 1 J 
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after some calculations wc get the symmetric hyperbolic system for the new vector of unknowns lA 
{q,u,h,S) (compare with (6), (44a)): 



/ 


Pp/p 





-{pp/p)h 





















('\ 






pao 


O3 


0^ 


dt 


u 


+ 


V 


03 


03 


0^ 




u 






O3 do - 


f {pp/p)h € 




h 




0^ 


03 


03 


0^ 




h 


V 











l) 








\o 








0^ 











(48) 



V 



u 
h 



( {PpIp)w-^ V- -{pp/p)hw-V 

V pcLQW ■ V — flo/i • V _ 

-{pp/p)h^w ■ V -do/i • V (do + {pp/p)h ® h)w ■ V 0^ 

w-Vj 

where do is the symmetric and positive definite matrix 

with a new matrix C in the zero-order term (whose precise form has no importance) and where wc have 
set F = di^i RF. We write system (48) in compact form as 

3 

AodtU + Y^{Aj + £ij+i)djU + C'U = F, 



(49) 



where 



£12 = 



( 


1 





••• 


^ 








( 





1 


• 


■ o\ 


1 








••• 




















• 


• 











••• 











1 








• 


• 











••• 


















• 


• 


V 








••• 


) 








I 








• 


■ 0; 










( 








1 ••• 


\ 



























••• 



















£14 = 











••• 



















1 








••• 





















I 








••• 


) 











The formulation (49) has the advantage of the form of the boundary matrix of the system Ai + £12 , with 

Ai = Q on WT, (50) 

because wi — hi = 0, and £12 a constant matrix. Thus system (49) is symmetric hyperbolic with 
characteristic boundary of constant multiplicity (see [16, 17, 18] for maximally dissipative boundary 
conditions). Thus, the final form of our reduced linearized problem is 

3 



AodtU + Y^{Aj + £ij+i)djU + C'U = F, 


in QJ, 


(51a) 


V X Jo = 0, div [) = 


in Q^, 


(51b) 


dt^p = ui - V2d2'p - vsdsip + (pdiVN, 




(51c) 


q = n-'H- [diolip, 




(51d) 


Un = d2{n2v) + d^in^ip) 


on lut, 


(51e) 


{u,n,ip) = o 


for t < 0, 


(51f) 



under the constraints (45), (46). 
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3. Function Spaces 
Now we introduce the main function spaces to be used in the foUowing. Let us denote 

Q±:=Rtxf]±, a;:=Rtxr. (52) 
3.1. Weighted Sobolev spaces. For 7 > 1 and s S M, we set 

and, in particular, A'*'^ :— . 

Throughout the paper, for real 7 > 1 and n > 2, i7^(R"') will denote the Sobolev space of order s, 
equipped with the 7— depending norm || • ||s^^ defined by 



\u 



|2 

Is, 7 



(2^)-" / A2^'^(e)|u(e)|2de, (53) 



u being the Fourier transform of u. The norms defined by (53), with different values of the parameter 7, 
are equivalent each other. For 7 = 1 we set for brevity || ■ ||s := || ■ (and, accordingly, the standard 
Sobolev space 7?*(K") := _ff|(IR")). For s S N, the norm in (53) turns to be equivalent, uniformly with 
respect to to the norm || • ||i/s(Ri) defined by 



ja I <s 

For functions defined over we will consider the weighted Sobolev spaces HJ^{Q^) equipped with the 
7— depending norm 

I Q I < m 

Similar weighted Sobolev spaces will be considered for functions defined on . 

3.2. Conormal Sobolev spaces. Let us introduce some classes of function spaces of Sobolev type, 
defined over the half-space Q^. For j = 0, . . . , 3, we set 

Zo = dt, Zi a{xi)di , Zj := dj , for j = 2, 3 , 

where ^{xi) € C°°(R-|-) is a monotone increasing function such that cr(xi) = xi in a neighborhood of 
the origin and (7{xi) = 1 for xi large enough. Then, for every multi-index a = (ao, . . . ,0-3) G N*, the 
conormal derivative Z" is defined by 

Z, .— ^0 • • • ^3 ' 

we also write d" = Oq" . . . d^'^ for the usual partial derivative corresponding to a. 

Given an integer ?n > 1, the conormal Sobolev space H"^^{Q^) is defined as the set of functions u € 
L'^iQr) such that Z"u e L^{Q^), for all multi-indices a with |a| < m (see [14, 15]). Agreeing with the 
notations set for the usual Sobolev spaces, for 7 > 1, Hl2n ^(Qt) ^^^^ denote the conormal space of order 
m equipped with the 7— depending norm 



i|<r 



and we have H12^{Q^) := i?t™„ liQr)- Similar conormal Sobolev spaces with 7-depending norms will be 
considered for functions defined on Q+. 

We will use the same notation for spaces of scalar and vector-valued functions. 
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4. The main result 

We are now in a position to state the main result of this paper. Reeah that lA = (g, u, h, S), where u 
and h were defined in (47). 

Theorem 8. Let T > 0. Let the basic state (22) satisfies assumptions (23) -(28) and 

\H y.'H\>5>Q onu^T, (55) 

where S is a fixed constant. There exists Jo > 1 such that for all 7 > 70 and for all J-~^ e -^tan fiQx)' 
vanishing in the past, namely for t < 0, problem (51) has a unique solution {Uj, Tij, (fij) G Htan ^(Qt) ^ 

HjiQx) ^ H}^{'^t) with trace {q^ , ui^ , hi^)\^,j^ G Hj^'^{ujt), H-yluiT ^ H-r^^i^r)- Moreover, the solution 
obeys the a priori estimate 

c 

where we have set U-^ = e^'^*U,'H-f = e^^*''H,Lp~^ = e""'* ip and so on. Here C = C{K,T,6) > is a 
constant independent of the data T and 7. 

The a priori estimate (56) improves the similar estimate firstly proved in [21]. 

Remark 9. Strictly speaking, the uniqueness of the solution to problem (51) follows from the a priori 
estimate (42) derived in [21], provided that our solution belongs to . We do not present here a formal 
proof of the existence of solutions with a higher degree of regularity (in particular, H^) and postpone this 
part to the future work on the nonlinear problem (see e.g. [4, 19]J. 

The remainder of the paper is organized as follows. In the next Section 5 we introduce a fully hyperbolic 
regularization of the coupled hyperbolic-elliptic system (51). In Section 6 we show an a priori estimate 
of solutions uniform in the small parameter e of regularization. In Section 7 wc show the wcU-poscdness 
of the hyperbolic regularization and in Section 8 we conclude the proof of Theorem 8 by passing to the 
limit as e — > 0. Sections 9, 10, 11 are devoted to the proof of some technical results. 



5. Hyperbolic regularization of the reduced problem 

The problem (51) is a nonstandard initial-boundary value problem for a coupled hyperbolic-elliptic 
system. For its resolution we introduce a "hyperbolic" regularization of the elliptic system (51b). We will 
prove the existence of solutions for such regularized problem by referring to the well-posedness theory for 
linear symmetric hyperbolic systems with characteristic boundary and maximally nonnegative boundary 
conditions [17, 18]. After showing suitable a priori estimate uniform in e, we will pass to the limit as 
e — >■ 0, to get the solution of (51). 

The regularization of problem (51) is inspired by a corresponding problem in relativistic MHD [22]. In 
our non-relativistic case the displacement current (1 / c)dtE \s neglected in the vacuum Maxwell equations, 
where c is the speed of light and E is the electric field. Now, in some sense, we restore this neglected 
term. Namely, we consider a "hyperbolic" regularization of the elliptic system (51b) by introducing a 
new auxiliary unknown which plays a role of the vacuum electric field, and the small parameter of 
regularization e is associated with the physical parameter 1/c. We also regularize the second boundary 
condition in (51d) and introduce two boundary conditions for the unknown E^ . 
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Let us denote = {H^,E'^). Given a small parameter e > 0, we consider the following regularized 
problem for the unknown {U^ ^V^ ,ip'^): 



3 

AodtW + Y^i^J + £i3+i)djU' + C'W = T 

£dt\f + V X £^ = 0, edtz^ - V x Jo^ = 
dt'V^ = It! - ^282^^ - V'id-iif' + ^'dxVN, 
q' =H- W - ~£E- E', 

E'^^_ ^edt{W)~ed2{E^ip'), 
E%^ = -edt{n2^n ~ed^{Ei^') 



in Qj,, 
in Q^, 



on wt, 
for t < 0, 



(57a) 

(57b) 
(57c) 

(57d) 
(57e) 
(57f) 
(57g) 



where 



E' ^iEt,EI,EI), E^iEi,E2,E3), r = (iJ^^i*!, i?;^, i?^, 
{E%,Eldi^i,Eldi^i), E%=El- Eld2^ - Eld3^, = E^dk^ 



El k 



-2,3, 



the coefficients Ej are given functions which will be chosen later on. All the other notations for "H"^ (e.g., 
S)"^ , f)^) are analogous to those for T-L. 

If ^I* = 0, $1 ^ xi, then f)" = J^*^ = T-L'^ , c'^ = = E^ , and when e = 1 (57b) turns out to be nothing 
else than the Maxwell equations. 

It is noteworthy that solutions to problem (57) satisfy 



+ 



div =0 in g 

div ()^ = 0, div e'' = in 

n% = d2{H2V>') + d3{n3^') on ujt, 



(58) 
(59) 
(60) 
(61) 



because (58)-(61) are just restrictions on the initial data which are automatically satisfied in view of 
(57g). Indeed, the derivation of (58) and (60) is absolutely the same as that of (45) and (46). Equations 
(59) trivially follow from (57b), (57g). Moreover, condition (61) is obtained by considering the first 
component of the first equation in (57b) at a;i = and taking into account (57e) - (57g). 

5.1. An equivalent formulation of (57). In the following analysis it is convenient to make use of a 
different formulation of the approximating problem (57), as far as the vacuum part is concerned. 

First we introduce the matrices which are coefficients of the space derivatives in (57b) (for e = 1 the 
matrices below are those for the vacuum Maxwell equations): 



Bl = £- 



( 





V 









-1 













1 




\ 

-1 





/ 



Bl = e~ 



( 








1 



1 
-10 

\ 



-1 









/ 









\ 





/ 







-1 






1 \ 









0/ 
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Then system (57b) can be written in terms of the "curved" unknown — (i^^, as 

3 

BodtW + B'jdjW + B4W' = 0, (62) 
i=i 

where 

Bo = idi$i)-^ KK^ >0, K^l2®f), Bi^dtBo, 
and the matrices Bq and K are found from the relations 

so that 

:.') - (t (s) - 

System (62) is symmetric hyperbolic. The convenience of the use of variables (i^', €^) rather than {W,E'^) 
stays mainly in that the matrices B^ of (62), containing the singular multiplier e~^, are constant. 

Finally, we write the boundary conditions (57c)-(57f) in terms of {W^,W'^), where we observe that 
(recalling that = 1 on ljt)'- 

n-w ^ 'Hn'hi + 112%%, + Tisn^ = ^ ^^ 

(63) 

E-E' ^ EnEI + E2E'^^ + i^gS;^, = ^ • €^ 

Concerning the first line above in (63) we notice that (ii='HAr = Oon cjt, so that does not appear 
in the boundary condition. 

From (62), (63) we get the new formulation of problem (57) for the unknowns {U^,W^): 

3 

AodtW + Y,(^3 + £i3+i)djU' + dW = J", in Q+, (64a) 

3 

Bodt W" + ^3 + B4W'' ^0 in Qj, , (64b) 

dtip" + V2d2ip'' + hdsv'' - f^diVN ~ul=0, (64c) 
+ [diq]ip'' + si-e^ ^0, (64d) 
€1 ~ e dtinsip^ + £ d2{Eiip') ^ 0, (64e) 

(Bl+edt{n2V')+ed3{Eiip') ouujt, (64f) 

(^^^ if') ^ for t < 0. (64g) 

From (58)-(61) we get that solutions {U',W^) to problem (64) satisfy 

dWh^=0 in Q+, (65) 

div ()^ = 0, div e'' = in , (66) 

hi = H2d2^' + H3d3^' - ^'diHN, (67) 

t)! = 82 {Ti2V') + 93 (H3</?") on UT. (68) 

Remark 10. The invertible part of the boundary matrix of a system allows to control the trace at the 
boundary of the so-called noncharacteristic component of the vector solution. Thus, with the system (64a) 
(whose boundary matrix is —£12, because of (50) J we have the control of q^,Ui at the boundary; therefore 
the components ofW^ appearing in the boundary conditions (64c), (64d) are well defined. 

The same holds true for (64b) where we can get the control ofSj2,^3, ^ii ^l- control of <t\ (which 
appears in (64d)j is not given from the system (64b), but from the constraint (66), as will be shown later 
on. We recall that ^\ does not appear in the boundary condition (64d) because \)i = T-Ln = 0. 
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Before studying problem (64) (or equivalently (57)), we should be sure that the number of boundary 
conditions is in agreement with the number of incoming characteristics for the hyperbolic systems (64). 
Since one of the four boundary conditions (64c)-(64f) is needed for determining the function ip'^(t,x'), 
the total number of "incoming" characteristics should be three. Let us check that this is true. 

Proposition 11. //O < e < 1 system (64a) has one incoming characteristic for the boundary lot of the 
domain Q^. If e > is sufficiently small, system (64b) has two incoming characteristics for the boundary 
ujt of the domain . 

Proof. Consider first system (64a). In view of (50), the boundary matrix on ujt is —£12 which has one 
negative (incoming in the domain Q^) and one positive eigenvalue, while all other eigenvalues are zero. 

Now consider system (64b). The boundary matrix Bf has eigenvalues A1.2 = — e^^, -^3^4 = e~^, A5.6 = 
0. Thus, system (64b) has indeed two incoming characteristics in the domain {^1,2 < 0). □ 

6. Basic a priori estimate for the hyperbolic regularized problem 

Our goal now is to justify rigorously the formal limit e ^ in (57)-(61), or alternatively in (64)-(68). 
To this end wc will prove the existence of solutions to problem (64)-(68) and a uniform in e a priori 
estimate. This work will be done in several steps. 

6.1. The boundary value problem. Assuming that all coefficients and data appearing in (64) are 
extended for all times to the whole real line, let us consider the boundary value problem (recall the 
definition of Q'^jU} in (52)) 



3 







in Q+, 


(69a) 


3 

BodtW + BpjW + BiW 


= 


in Q", 


(69b) 


dtip" + V2d2^^ + V'ids.tp^ - V^diVN ~ u\ - 


= 0, 




(69c) 




= 0, 




(69d) 


^l~edt{nzip')+ed2{Ei^') : 


= 0, 




(69e) 


+£at(H2¥'") +e 53(^1 (^") 


= 


on oj, 


(69f) 




= 


for t < 0. 


(69g) 



In this section we prove a uniform in e a priori estimate of smooth solutions of (69). 

Theorem 12. Let the basic state (22) satisfies assumptions (23)~(28) and (55) for all times. There exist 
£0 > 0, 70 > 1 such that if < e < eo and 7 > 70 then all sufficiently smooth solutions {U^ ^Lp^) of 
problem (69) obey the estimate 

+ 7^||v:';il?,.(.)<^||-F,|l?,._(^,), (70) 

where we have set lA^ = e^'^^U^ , = e~'^^ W' , (ys^ = e~'^* (p' and so on, and where C — C{K, S) > is 
a constant independent of the data J- and the parameters £,7. 

Passing to the limit £ — >■ in this estimate will give the a priori estimate (56). 

Since problem (69) looks similar to a corresponding one in relativistic MHD [22], for the deduction 
of estimate (70) we use the same ideas as in [22]. On the one hand, we even have an advantage, in 
comparison with the problem in [22] , because the coefficients Ej in (69b) , (69d)-(69f) are still arbitrary 
functions whose choice will be crucial to make boundary conditions dissipative. On the other hand, we 
should be more careful with lower-order terms than in [22], because we must avoid the appearance of 
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terms with e ^ (otherwise, our estimate will not be uniform in e). Also for this reason we are using the 
variables {U^,^) rather than (U^^V). 

For the proof of (70) we will need a secondary symmetrization of the transformed Maxwell equations 
in vacuum (57b), (59). 

6.2. A secondary symmetrization. In order to show how to get the secondary symmetrization, for 
the sake of simplicity we consider first a planar unperturbed interface, i.e. the case ip = 0. For this case 
(57b), (59) become 



(71) 



div-H'=0, div£;^ = 0. (72) 
We write for system (71) the following secondary symmetrization (for a similar secondary symmetrization 
of the Maxwell equations in vacuum see [22]): 

3 

%dtV' + J2 %B'jdjV' + RiAwW + i?2divS" 



0, 



(73) 



where 



/ 1 





£1^3 






1 









V 




£1/2 




£1/1 

















































-1^1 























V\ 




V2 




v-i 








e-i 






-v\ 









V 
















-V, ) 






/ ~vz 


V\ 







-£-1 \ 











V2 


£ 


-1 








%% = 





£-1 


V-i 













1^1 




—e 


-1 










-J/g 


V2 




\ 










V2 


1^3 J 



£1/3 — £1/2 \ 
~£I/3 £i/i 
£1/2 — £J/l 

1 
1 
1 / 

/ -1/2 J/1 
J/1 J/2 1/3 
J/3 -J/2 

-£-- 


V £-1 

/ ^1 \ 

J^2 
J^3 




V J 



(74) 



Ri = 






-J^2 
J^l 





R2 = 



£-1 \ 


j/i 

J^2 J^3 

J/3 -1/2 y 

/ \ 




j^i 

J/2 
V ^3 / 



The arbitrary functions Vi(t,x) will be chosen in appropriate way later on. It may be useful to notice 
that system (73) can also be written as 



{dtW + X E") - i7 X {edtE" - V x H^) + vdivW 

£ 



0, 



1 



(75) 



V X ■H'') + j7x {edtW + V x E') + iydivE^ = 0, 



with the vector-function z7 — (j/i, j/2, J^s). The symmetric system (73) (or (75)) is hyperbolic if *Bq > 0, 
i.e. for 

£|j7| < 1. (76) 
The last inequality is satisfied for any given j/ and small e. We compute^ 

detCBD^^iMl'^l'-lA')'- 



-'^The manual computation of the determinants is definitely too long. Here we used a free program for symbolic calculus, 
with the help of PS's son Martino. 
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Therefore the boundary is noncharacteristic for system (73) (or (75)) provided (76) and i^i 7^ hold. 

Consider now a nonplanar unperturbed interface, i.e., the general case when Lp is not identically zero. 
Similarly to (73), from (62), (59) we get the secondary symmetrization 



K'^IK-^ \ B^dtW^ + B^^d.W^ + BiW^\ + ^^A'fi?idiv fj^ + i?2div ( 



= 0. 



We write this system as 



where 



3 

M^dtW^- + M^dj W + MIW" =Q, (77) 

= -L- KWr.K'^ > 0, An = KW.K'^ (j = 2, 3), 

1 ~ ~ 1 ^ 

Mf^-^K^IK^, ^1= f^l-Y^^ld,^)^ (78) 



System (77) is symmetric hyperbolic provided that (76) holds. We compute 

det(A/f) = (1 + (52^)2 + (33^)2)' {lyi ~ ^282^ - i^sd^ipf - l/e'f , (79) 

and so the boundary is noncharacteristic for system (77) if and only if (76) holds and vi 7^ iy2d2'P + V3d3ip. 
System (77) originates from a linear combination of equations (57b) similar to (75), namely from 



(dti)" + - V X €^) - ?7(i7x fj-^edte" -V x S)^)) + ^^divf)'' = 0, 
£ i9i$i 

{dtz' - iv X Sj^) + fj (v X 7rHe<9tf)^ + V x r)) + dive^ = 0. 



(80) 



We need to know which is the behavior of the above matrices in (78) w.r.t. e as e — >■ 0. In view of this, 
let us denote a generic matrix which is bounded w.r.t. e by 0(1). Looking at (80) we immediately find 

K = 0(1), MJ = 5^ + 0(1) (i = 1,2,3), Ml^O{l). (81) 

As the matrices Mg and do not contain the multiplier £~^, their norms are bounded as e — >■ 0. 
Recalling that the matrices i3| are constant, we deduce as well that all the possible derivatives (with 
respect to t and Xj) of the matrices have bounded norms as £ — 0. 

6.3. Proof of Theorem 12. For the proof of our basic a priori estimate (70) we will apply the energy 
method to the symmetric hyperbolic systems (69a) and (77). In the sequel 70 > 1 denotes a generic 
constant sufficiently large which may increase from formula to formula, and C is a generic constant that 
may change from line to line. 

First of all we provide some preparatory estimates. In particular, to estimate the weighted conormal 
derivative Zi = crdi of (recall the definition (54) of the 7-dependent norm of H^^^^ ^) we do not need 
any boundary condition because the weight a vanishes on w. Applying to system (69a) the operator Zi 
and using standard arguments of the energy method^, yields the inequality 

7ll^imii.(Q.) < ^ {m\lu.,iQ^) + ll"7llk„,.(Q+) + ll«7lli^(Q-)} ' (82) 
for 7 > 7o. On the other hand, directly from the equation (69a) we have 

\\£i2u%.^Q,^ < c{ii + mr„i^^_^^Q.^} , (83) 



^ We multiply Zi(69a) by e Z\US. and integrate by parts over Q+, then we use the Cauchy-Schwarz inequality. 
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where C is independent of £,7. Thus from (82), (83) we get 

7ll^imii^(Q+) < ^ + ll^7llk„,.(Q+)} ' ^ ^ (84) 

where C is independent of e, 7. Furthermore, using the special structure of the boundary matrix in (69a) 
(see (50)) and the divergence constraint (65), we may estimate the normal derivative of the noncharac- 
teristic part U^^ — e~^^{(f , Ui,hl) of the "plasma" unknown U^: 

}, (85) 

where C is independent of e, 7. In a similar way we wish to express the normal derivative of through 
its tangential derivatives. Here it is convenient to use system (69b) rather than (77). We multiply (69b) by 
e and find from the obtained equation an explicit expression for the normal derivatives of ^331 -^l' 'SI' ^1- 
An explicit expression for the normal derivatives of iof , €f is found through the divergence constraints 
(66). Thus wc can estimate the normal derivatives of all the components of through its tangential 
derivatives: 

WdiW^WliQ-) < C |7^I|W^,1IL(Q-) + WdtWX.^^,) + E , (86) 

where C does not depend on e and 7, for all £ < Eq- 

As for the front function (p'^ we easily obtain from (69c) the estimate 

7 

where C is independent of 7. Furthermore, thanks to our basic assumption (55)'^ we can resolve (67), 
(68) and (69c) for the space-time gradient Vt,a;'<P^ = (9ti^^, 92</?^, c^s'/'^): 

^t,x'H^% = axh\^ + a2f)i^ + azu\^ + a^ipl^ + ^a^tpl^, (88) 

where the vector-functions = aa{U\^^'H\^) of coefficients can be easily written in explicit form. From 
(88) we get 

I|Vm'(^^IU.(^) < c {Wu^^Wl^^) + Ww^IWl-^^) + 7I1v^IU^m) ■ (89) 

Now we prove a energy estimate for (U^ , W^). We multiply (69a) by e~'^^U^ and (77) by e"''* W^, 
integrate by parts over Q*. then we use the Cauchy-Schwarz inequality. We easily obtain 



iMJh^^^ < -WulfL^u), 7>7o, (87) 



' / {AoU^,U^)dxdt + -f [ {M^W^,W^)dxdt+ [ A' dx'dt 

Jq+ JQ- Ju) 



e\\2 



(90) 



where we have denoted 

1 1 



A' = --(£l2^,^)U + o(^-^l'^7''^7)l^ 



2 

Thanks to the properties of the matrices {a = 0,4) described in (81), the constant C in (90) is 
uniformly bounded in e and 7. Let us calculate the quadratic form A"^ for the following choice of the 
functions Vj in the secondary symmetrization"*: 

vi=V2d2'p + v3d3(p, Vk = vu, fc = 2,3. (91) 

After long calculations we get (for simplicity we drop the index 7) 

= -g^?/f +£-l(-^3^2"-*l2*S3) + («2.^|+«3-^|)H^ + (t)2£|+V3C|)£^^, On W. (92) 



■^Under the conditions //jv = 'Hjv = one lias \H X "Hp = {H2'H3. - -f^3'H2)^ ( V^)^ on uj, wiierc we have set (V'l^) : = 
{l + \d2^\^ + \a:i^??/-'. 

^Notice that the choice (91) makes the boundary characteristic, see (79). 
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Now we insert the boundary conditions (68), (69c)-(69f) in the quadratic form A'^ , recaUing also HnIu = 
and noticing that 

Again after long calculations we get 

A' = {El + «2H3 - hn2) {sE'^dt^^ + S^ld^iip' ~ ^ld2v') 

+ {eEr.El + eEr.El) {dtv' + V2d2v' + hd^v") 

+ - jq' + [diq] ul - diVNiq' + [diq]^') + {lU^ + OM^ - d2Ei){Sjl + ehEf,) 

+ {in2 + dtn2 + d:iEi){^ll~ev:iEI,) + {d2n2 + d:inz){v2m + ^z^l)) on cj . (93) 

Thanks to the multiplicative factor e in the boundary condition (69e), (69f), the critical term with the 
multiplier e^^ in (92) has been dropped out. We make the following choice of the coefficients Ej in the 
boundary conditions (69d)-(69f): 

E^ilxiy, 

where ? is that of (91). For this choice 

El + V2n3 - = 0, Er, = 0, Er, = 0, (94) 

and this leaves us with 

A"" = (p-{ - -f{q^ + [diq](p-) + [diq] {ul + Lp^divn) + diVNq" 

+ (7H3 + dtU^ - d2Ei){^l + ev2El,) + (7H2 + dtHi + 93^i)(-^2 - ^^El,) 

+ (92^2 + 93^3 )(l'2i32 + l'3i03)} LO . 

Then we write in more convenient form the terms with coefficient 7 substituting from (69d) 

-(q^ + [diq]^') + n2S)i + = £ e • r , 

and we notice that 

e • r + {V2'H?. - V3'H2)EIj = e • + (^2^3 - ^112^1 = 0, on w, 
again by (94). Thus we get 

A" = ^"{[diq] {u\ + ip^diVN) + diVNq" 

+ {dtUz " d2Ei){^f^ + ev2E%) + {8^2 + ^3^1) (-^2 - e^^E^) 

+ {d2n2 + d3n3){v2Sjl + h^l)} on uj. (95) 

From (90), (95) we obtain (we restore the index 7) 

+ C (lK1li.(Q.) + \\W%.^Q-)) + l\m\hi.)^ (96) 
where C is independent of £,7. Thus if 70 is large enough we obtain from (87), (96) the inequality 

^ f {ll-^7!li^(Q+) + WK-rUhi^) + WW^UIh.)}^ < 6 < £0, 7 > 70, (97) 

where C is independent of £, 7. 
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Now we derive the a priori estimate of tangential derivatives. Differentiating systems (69a) and (77) 
witfi respect to Xq = t, X2 or xs, using standard arguments of tlie energy method, and applying (85), 
(86), gives the energy inequality 



7 / {AoZeU^,ZiU;)dxdt + -f / {M^ZeW^, ZeW^)dxdt + / Afdx'dt 



7 

where £ ~ 0,2,3, and where we have denoted 

1 



A} 



-(£i2Zim,Zem 



-(MfZeW'ZiW^)\^. 



m\\h,iQ-)}^ (98) 



Thanks to the properties of the matrices {a = 0,4) described in (81), the constant C in (98) is 
uniformly bounded in e and 7. We repeat for Al the calculations leading to (95) for A"^ . Clearly, for the 
same choices as in (91) and (94) we obtain (for simplicity we drop again the index 7) 

Al = Zeif^lldiq] {Ziul + Znp^diVN) + diVNZiq"^ 

+ {dtUz - d2Ei){Zi<ol + shZiEff) + {dt'H2 + d:iEi){Z,^l - ev^ZiEf,) 

+ (aaTis + d3n3){v2ZiSjl + V3ZiF)l)} + l.o.t., on uj, (99) 

where l.o.t. is the sum of lower-order terms. Using (88) we reduce the above terms to those like 

chlZiul, ch\Ziip^, chlZiS^^j, c.h\Zi(t'^j, ... on w, 

terms as above with f)f,uf instead of /if, or even "better" terms like 

-fCip^Ziu\, ^Cip'^Z^Lp" . 

Here and below c is the common notation for a generic coefficient depending on the basic state (22). By 
integration by parts such "better" terms can be reduced to the above ones and terms of lower order. 
The terms like ch\Z(U^^^^Q are estimated by passing to the volume integral and integrating by parts: 

/ ch\Ziu\,^^^QAx' At = — I di{ch\Ztu\)(ixAt 
Juj Jo+ 

^{Zid)hl{diul) + d{Zihl)diul - {did)hlZiul-c{dihl)Zeuljdxdt, 

where cj^^j^o = c. Estimating the right-hand side by the Holder's inequality and (85) gives 

/ chlZ,ul^^^^,dx'dt <c{||^^||i.(Q+) + |m||^,^^ (100) 

In the same way we estimate the other similar terms c h\ZiSj'^j , c h\Zi(B'!j , etc. Notice that we only need 
to estimate normal derivatives either of components of U^^ or W^. For terms like c f)f Z^uf , c f)f etc. 
we use (86) instead of (85). 

We treat the terms like c/if j^^^gZ^cp'^ by substituting (88) again: 



/ ch\Zg(p'^ dx' dt 

J U) 



ch\{dih\ + d2i)l + a^ul + d^ip^ -f 'yd^^ip'^ \dx dt 



< ciimuh^^^ + liw^1c.|li.(.) +7ll^1li2M . (101) 



PLASMA- VACUUM INTERFACE 



21 



Combining (98), (100), (101) and similar inequalities for the other terms of (99) yields (we restore the 
index 7) 

7 (ll^^mii^tQ.) + WZiW^.^Q-)) < c^{;ll-^7llk„,.(Q+) + ll"7llk„,.(Q+) + W^^WkiQ-) 

+ 7 {\\K^Uh(^^ + WW^UU^)) }, < £ < £0, 7 > 70, (102) 
where C is independent of £,7. Then from (84), (86), (97), (102) we obtain 

+ l{\\KM\hi^) + \\W^U\h^u))}, 0<£<£o, 7>7o, (103) 
where C is independent of e, 7. We need the following estimates for the trace of U^, W^. 
Lemma 13. The functions Uf^jW^ satisfy 

7l|Z^,t,Ulli.(„) + m.Uly^^^^ < C (ll-^.lli.(Q.) + ll"7llk„,.(Q+)) ' (104) 

l\\W^U\l.^^) + mUlu^^^^ < C\m\\h^^Q-y (105) 

The proof of Lemma 13 is given in Section 11 at the end of this article. Substituting (104), (105) in 
(103) and taking 70 large enough yields 

^ + W^^WkiQ-)) ^ 7ll-^^llk„,.(Q+)' < e < £0, 7 > 70, (106) 

where C is independent of £,7. Finally, from (89), (104) and (106) we get 

Adding (106), (107) gives (70), The proof of Theorem 12 is complete. 

7. WeLL-POSEDNESS of the hyperbolic REGULARIZED PROBLEM 

In this section we prove the existence of the solution of (69). Its restriction to the time interval (—00, T] 
will provide the solution of problem (64). From now on, in the proof of the existence of the solution, e is 
fixed and so we omit it and we simply write U instead ofW, W instead of W^, (p instead of ip^ . 

In view of the result of Lemma 16 (see Section 9) we can consider system (77) instead of (69b). First of 
all, we write the boundary conditions in different form, by eliminating the derivatives of tp. We substitute 
(69c) in the boundary conditions for £2, (£3 and take account of the constraint (68) and the choices (91), 
(94). We get 

q - f)2i32 - ^3-^13 + e EiEn + [diq]Lp = 0, 

€2- eilaui+evaHN +eaiip ^0, (108) 
£3 + £"^21*1 — £'C2'Hjv + ea2iy9 = 0, on w, 

where the precise form of the coefficients ai,a2 is not important. For later use we observe that (68), 
(69c)-(69f) is equivalent to (68), (69c), (108). Notice that the last two equations in (108) yield 

e Eiui + ^2^2 + vs'Bs + eastp = 0, (109) 

where 03 = aiV2 + 02^3. 
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Let US write the system (69a), (77), (108) in compact form as 

C 



U 
W 
U 

w ^ 
t(W,W^,^) = 







M 



+ 6 </? = 0, in 

for t < 0, 



(110) 



where the matrix M and the vector b are imphcitly defined by (108). 

Let us multiply (110) by e"'''* with 7 > 1; according to the rule e~'^*dtu — {j + dt)e~'''^u, (110) becomes 
equivalent to 



M 



+ b(p~f =0 



on (3+ X 

in cj, 
for t < 0. 



(Ill) 



where 



^0 

Aq 



e~^^U, = e-''* W, = e-^V, etc. 



First we solve (111) under the assumption that </?^ is given. 

Lemma 14. There exists 70 > such that for all 7 > 70 and for all given T G e"'*H^^^ ^{Q^) and 
tp G eJ*Hj^^{Lu) vanishing in the past, the problem (111) has a unique solution {U, W) G e'^*Hlg^„ ,y{Q^) x 
e'^^HUQ") with {q,uihi,W)\^ G 6^*77^^(0;), such that 



|le-^*(Z^,M^)|Ui__^(Q+),^.(Q-) + i|e-^*(g,7.i,/ii,M^)|„|l^v.(^^ 



7 



(112) 



Proof. We insert the new boundary conditions (108), (109) in the quadratic form (see (92)) and we 
get 

A' -.^ -]^{Ai+£i2)U- U + ]^MIW -W ={[dq]ui +02^2- ai^z-ea:,Ej^)^ oncj. (113) 

If we consider the boundary conditions (108), (109) in homogeneous form, namely if we set tp = 0, then 
from (113) 

= on uj. 

We deduce that the boundary conditions (108) are nonnegativc for C-y. As the number of boundary 
conditions in (108) is in agreement with the number of incoming characteristics for the operator C-y (see 
Proposition 11) we infer that the boundary conditions (108) are maximally nonnegative (but not strictly 
dissipative) . Then we reduce the problem to one with homogeneous boundary conditions by subtracting 
from (U-y.W-y) a hmction (Z^^, W!^) G H^iQ+) x H^iQ') such that 



M 



W 
W 



b(p ^ on oj. 



Finally, as the boundary is characteristic of constant multiplicity [16], we may apply the result of [17, 18] 
and we get the solution with the prescribed regularity. □ 



The well-posedness of (69) in H^^^ x is given by the following theorem. 
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Theorem 15. There exists 70 > such that for all 7 > 70 and J- G e^^H^^^^ ^{Q'^) vanishing in the 
past, the problem (69) has a unique solution {U, W) G e^^Hl^^^ -yiQ^) ^ ^^^H}^{Q~) with {q,uihi, W)\^ G 

Proof. We prove the existence of the solution to (69) by a fixed point argument. Let Ip g e'^*H^^'^{u!T) 
vanishing in the past. By Lemma 14, for 7 sufficiently large there exists a unique solution (U, W) € 
eT*i/i,„,^(g+) X 6-^*771 (g-), with {q,ui,hi,W)i^ e eT*i?^/'(a;) of 



M 



on Q+ X Q , 



(114) 



_(iY^,VF^) = fort<0, 
enjoying the a priori estimate (112) with Ip instead of (/S. Now consider the equation 

7V3^ + 9f93^ + -0232<P7 + — <y579iujv = tti7, onw, (115) 

1 /2 

where U17 G ^7 (uS) is the trace of the component of given in the previous step, vanishing for t < 0. 

1/2 

For 7 sufficiently large there exists a unique solution 1^7 G H~f (cj), vanishing in the past, such that 

C, 



\Wi\\hI 



< —\\u 
7 



(116) 



From the plasma equation in (114) and from (115) we deduce the boundary constraint 

/117 — -ff2(?2</?7 + H^d^ipj — ipjdiHN on w. (117) 

Since in the right-hand side of (114) we have Ip instead of (p we are not able to deduce the similar boundary 
constraint for the vacuum magnetic field. Instead, we obtain 



f)i7 - d2{T-L2fp^) - 93(7^3(^7) = G7 on 



(118) 



where G-y solves 



—07 + a292(',f7 - V7) ~ ai93((p7 - ip^^) + {8202 - dzai){ip^ - ip^) = 



on w. 



(119) 



for d/dt = 7 + 9t + d2{v2-) + d3{v3-) and where the coefficients ai,02 are the same of (108). (119) is 
derived from the first equation of the vacuum part of (114), (115) and the boundary conditions for €2, £3 
in (114). 

Let us consider the linear system for Vt.x'^i provided by equations (115), (117) and (118). By the 
stabffity condition (55) we can express \7t,x''fi'y through {hij,t)i^,uij)^^, ip^ 



^t,x"f'i = a[hij + a2f)i7 + 031417 + a'^ipj + a'^Gj, 
where the precise form of the coefficients has no interest. Then, substituting into (119) yields 



— G7 + 60G7 = bihij + 62f)i7 + bsip^ + a2d2p^ - aidsp^ + {8202 - d^aijip^ 

with suitable coefficients 6,. 

From (121), for 7 sufficiently large, we get the estimate 

G / 



on W, 



< 



c 



{\\^i\\hI 



(120) 
(121) 



(122) 
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where we have apphed (112) (with Lp in place of ip) and (116). Thus, from (120) again, we obtain the 
estimate 



c 

< 

1 



C / \ 



Combining (112) (with p in place of Lp), (116) and (123) gives 

C_ 

7 

This defines a map ^ ^ ^ in e^^H^''^{uJT)- Let^\^^ e e^^ H^'"^ {lot), and {U^,W^), {U^,W^), V^,^^ be 
the corresponding solutions of (114), (115), respectively. Thanks to the linearity of the problems (114), 
(115) we obtain, as for (124), 

C 

Then there exists 70 > such that for all 7 > 70 the map Tp ^ has a unique fixed point, by the 
contraction mapping principle. The fixed point Tp = ip, together with the corresponding solution of (114), 
provides the solution of (111), (115), that is a solution of (69). As for the boundary conditions, we have 
already observed that (68), (69c)-(69f) is equivalent to (68), (69c), (108). The proof is complete. □ 

8. Proof of Theorem 8 

For all e sufficiently small, problem (64) admits a unique solution with the regularity described in 
Theorem 15. Due to the uniform a priori estimate (70) we can estract a subsequence weakly convergent 
to functions {lA, W, ip) with (Uj,Wj) G HI^^ ,^{Q1^) x H}^{Q:^) and (97,^1^, hi^)\^^ £ hI^^'^ {ut), W^\u,t € 
H}^{(jJt) and ^p-y £ H}^{ojt) (we use obvious notations). Let us decompose W — (.$3,(2) and perform a 
inverse change of unknown with respect to that of Section 5.1 to define ("H, E) from (io, £). Passing to 
the limit in (57b), (64)-(68) as e — )■ immediately gives that {11,%, ip) is a solution to (51), (45), (46) 
and E ^ >B — Q. Passing to the limit in (70) gives the a priori estimate (56). The proof of Theorem 8 is 
complete. 

9. Equivalence of systems (57b) and (80) 

We prove the equivalence of systems (57b) and (80) for every v ^Q. This is the same as the equivalence 
of (64b) and (77), or (69b) and (77). 

Lemma 16. Assume that systems (57b) and (80) have common initial data satisfying the constraints 

div \]^ = 0, div = in for t = 0. 

Assuming that the corresponding Cauchy problems for (57b) and (80) have a unique classical solution on 
a time interval [0,T], then these solutions coincide on [0,T] for all e sufficiently small. 

Proof. Let us set 

A = fi~\dt\)' + s-^V X r), B = fi-^dtc" - e^^V x Sj^). 
Then (80) can be written as 

A-ei^x B + ^^divi)' = 0, B + e v x A+ div e"" = 0. (125) 

Taking the vector product of t7 with the systems in (125) gives 

i7xA-£i7x(i7xB)=0, iyxB + eh'x{iyxA)=0, (126) 

that is 

i7x A- £{i7- B)i7 + e\iJ\'^B ^0, x B + e {i7 ■ A)i7 - e\i7\'^A ^ 0. (127) 
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We take the vector product oi ev with the first system in (127) and get 

e{v- A)v~e \v\^A + \v\^ z7 x B = 0. (128) 

For any choice of i7 7^ we may assume that e\v\ 7^ 1 (this is true for e definitely smaU) . Then by 
comparison of (128) and the second equation in (127) we infer v x B ~ and from (126) also v x A~ Q. 
Thus (80) may be rewritten as 

dt^" + iv X r + div = 0, dtt" xS^" + div = 0. 

Applying the div operator to the equations gives the transport equation 

dfU + div(wa) = in Q^, 

for both u = div f)^ and u = div e"^, where a = fjiJ/di^i. Noticing that the first component of a vanishes 
at Xi = 0, the transport equation doesn't need any boundary condition. As M|t=o = 0, by a standard 
argument we deduce m = for t > 0. This fact shows the equivalence of (57b) and (80). 

□ 



10. Proof of Lemma 1 

Given an even function x G Co° (^) ' with x = 1 on [—1,1], we define 

^{xi,x'):=xixi{D))ip{x'), (129) 

where x{xi{D)) is the pseudo-difFcrential operator with (D) = (1 + |Z?p)^/^ being the Fourier multiplier 
in the variables x' . From the definition it readily follows that ^'(0, x') ~ f{x') for all x' e R^. Moreover, 

d^-^{xi,x')=x!{xi{D)){D)^{x'), (130) 

which vanishes if cci = 0. We compute 

ii*(^i,-)iIh™(m^)- / {i'?"'x\x,{i'))m')?d£,' , 

where <,5(^') denotes the Fourier transform in x' of ip. It follows that 

In a similar way, from (130), we obtain 

Iterating the same argument yields 

ll^i*lli2^(H'"-j(R2)) < C|l(^||^„-o.5(R2) , j = 0,...,m. 



Adding over j = 0, . . . ,m finally gives ^ G H"^{W) and the continuity of the map ip 1-^ 

We now show that the cut-off function x, and accordingly the map Vf, can be chosen to give (14). 
From (130) we have 

d^^{x,,x')^{2n)-' [ e^«'-'x'(^i(e'))(CV(C')de'- 
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By the Cauchy-Schwarz inequality and a change of variables we get 

1/2 



We change variables again in the integral above by setting s = xi{p). It follows that 

\d^^{x)\<cMH^R^)[^J^ \x'i^)\'^-) <cMHHm[j^ '^'("^I'tJ ■ ^^^^^ 

Given any M > 1, we choose x such that x(s) = for \s\ > AI, and |x'(s)| !i 2/A/ for every s. Then 
from (131) one gets 

(J 



Given any e > 0, if M is such that C/^/M < e, then (14) immediately follows. 

The proof of Lemma 2 follows from Lemma 1, with t as a parameter. Notice also that the map (/3 — > ^P, 
defined by (129), is linear and that the time regularity is conserved because, with obvious notation, 
^'(9j(p) = df 'i>{Lp). The conclusions of Lemma 2 follow directly. 



11. Proof of Lemma 13 



We write U^^^ on lj as 



which gives 



/>oo 



WK^Uh^^^ < 2|1Z^^-|U.(Q+)||W,^|U.(Q+). (132) 
Estimating the right-hand side of (132) with (85) and using the 7-liomogeneity of the H^^^ ^ norm gives 

Thus the first part of (104) is proved. To show the second part of (104) we compute for i ~ 0,2, 3, 
/ \ZeU'„^fU,=o dx'dt = -2 / / ZiUl^ ■ diZdAl^ dxdt = 2 / / Z^U'„^ ■ diU^^ dxdt, 

J LO Jo J LO Jo J UJ 

which gives 

\m^U\]J^^^^) < 2\\U^\\H2^^^^^Q+^\\diUmL^Q^y (133) 
Interpolating between (132) and (133) gives 

< ni^',\\HU^iQ + )\\dlK,\\LHQ+)- 

Applying (85) eventually gives the second part of (104). We do the same for (105). 
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